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Abstract: The tau hadronic width provides a determination of the strong coupling con- 
stant a s at low energies, since it can be related to a weighted integral of the Adler function 
in the complex energy plane. Using Operator Product Expansion, one sees that the sensitiv- 
ity to a s comes from the perturbative contribution, which can be obtained by integrating 
the perturbative expansion of the Adler function. Two different prescriptions proposed 
to perform this integral, called Fixed-Order Perturbation Theory and Contour-Improved 
Perturbation Theory (FOPT and CIPT), yield different results for the strong coupling 
constant. Recently, models for the Adler function based on renormalon calculus have been 
proposed to determine which of the two methods is the most accurate, by comparing the 
resulting asymptotic series with the true value of the integral. We discuss the assumptions 
of such ansatz and the determination of their free parameters. We show that variations of 
this renormalon ansatz can yield opposite conclusions concerning the comparison of CIPT 
versus FOPT, and that such models are not constrained enough to provide a definite answer 
on this issue or to be exploited for a high-precision determination of a s {m^). 
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1. Introduction 

A precise assessment of the fundamental parameters of the Standard Model is mandatory 
to test the consistency of the theory and to exploit its full predictive potential, with the 
aim of identifying discrepancies indicating loopholes in our understanding of its dynamics 
or providing signs of New Physics. The coupling constant a s is the central ingredient of 
the strong sector of the theory, which can be determined at different scales using a large 
range of processes At low energies, a precise value of as(M^) can be extracted from 
the r hadronic width R T 

= gjr ->• hadrons v T (j)] 

T[t -»• eu e v T (j) K ■ ' 

as well as its moments [Q, [|, which can be determined from the analysis performed by the 
LEP experiments Q, and complemented with further experimental input from £?-factories. 

The theoretical description of this quantity can be obtained by relating this decay 
width with an integral of the r spectral functions, corresponding to the imaginary part 
of two-point correlators of hadronic currents. Exploiting the analytic properties of theses 
correlators, one can reexpress the integral as a contour integral of the Adler function D 
along a circle of radius |s| = in the complex plane of center-of-mass energy. One can 
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then use the Operator Product Expansion of the correlator || to compute the various 
contributions to the decay width |2], ^ || [7], || : 

Ky/A = ^S EW \V ud \ 2 [l + +6' EW +J2 C]v/a\ (1-2) 

D>2 

with Sew collects short-distance electroweak corrections || [[(]] , 5^ stems from the purely 
perturbative part of the expansion, defined in the chiral limit, while mass corrections are 
incorporated in (higher) D-dimensional contributions y/J containing also the contri- 
butions from the condensates. 



The main sensitivity to a s comes from 5^ , defined as 

-w(s)D(s), w(s) = 1-2- + 2 (- ,3 

\s\=sq 



1 + 5 <P) = _ 27Ti I ^- w ( s )D(s), = 1 - 2— + 2 f —V - f-V (1.3) 

J\s\=s S S \S J \S J 



Since we know only the first orders of the perturbative expansion of the Adler function for 
a given (real) value of the scale, the method used to compute the contour integral turns 
out to be important, in order to control (and suppress) the size of unknown higher orders. 
In particular, two different rules have been proposed, called Fixed-Order Perturbation 
Theory (FOPT) and Contour-Improved Perturbation Theory (CIPT) |], |7[ ||. We will 
recall below the most salient features of the two approaches. It turns out that the increasing 
experimental accuracy and the determination of the 0{cQ) term in the perturbative series 
of the Adler function has proved the two different approaches to differ significantly and 
to induce a significant systematic uncertainty in the determination of the strong coupling 
constant. 

Several studies have been performed to determine which of the two methods (if any) 
is to be preferred. A first ("internal") way of dealing with this issue consists in comparing 
the two methods to determine if one has a more regular and stable behaviour than the 



other, proving thus its robustness. For instance, in ref. [11|, our current knowledge of the 
perturbative series for the Adler function led us to conclude that CIPT showed a much 
better stability than FOPT as far as the dependence of the scale defining the integration 
contour is concerned, and that the FOPT integrand showed a pathological behaviour once 
one gets close to the end of the integration circle s = M%. In ref. [12], a similar study was 



performed, where the reference point for the Taylor expansion of the coupling constant used 
in FOPT was varied on a circle of radius M T in the complex energy plane, hinting at a 
strong dependence of the FOPT value due to a large impact of logarithmic corrections in the 
Taylor expansion. However, let us stress that these conclusions rely on the determination 
of a s at a reference point in the complex plane by iterating the RGE step by step along 
the contour of integration (in a similar way to the CIPT method), and thus use a "non- 
canonical" version of FOPT including elements of the CIPT philosophy. A model for higher 
orders in the perturbative expansion of the Adler function suggested that FOPT could 
actually oscillate towards the (more stable) CIPT value once higher-order are included. 

One can also opt for a different ("external") approach, where the true value of the 
Adler function and its perturbative expansion are assumed to be known. One can then 
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determine if the true value of the integral is approached by one of the two methods when 
one starts increasing the accuracy of the perturbative expansion (before the series becomes 
asymptotic). A particular model was proposed more recently in ref. Jl^| , based on renor- 



malon calculus |l3|, [T(], [17|, 15]. The first observation consists in the fact that perturbative 
series like 

oo 

D(Q 2 ) = 1 + J2 c ^ a Q a Q = as(Q)/ir (1-4) 

n=l 

are at best asymptotic ones, with a zero convergence radius in oq, but their Borel transform 
can be defined as 

B{B](t)=J2cn+i,i^ Tri (1.5) 
z — «■ 7r n+i n! 

n=0 

with improved convergence properties, and in particular a non-vanishing convergence ra- 
dius. If B[D](t) has no singularity for t real and positive, and does not increase too quickly 
at infinity, one can define the Borel sum: 

/•oo 

D(a) = dt e*/ a 5[D](t) (1.6) 
Jo 

The D(a) has the same perturbative expansion in powers of a as D(Q 2 ) in a s (Q). Actually, 
the Borel transform is expected to have singularities along the real axis, for borh positive 
and negative values of t. The former are called infrared renormalons, and are related 
to power corrections and condensates in the Operator Product Expansion of the Adler 
function, whereas the latter are called ultraviolet renormalons and are related to the large- 
order (and often oscillatory) behaviour of the series at higher orders. The presence of 
infrared renormalons requires to give a prescription to avoid the singularities in eq.|1.6|) 
(most often the principal value). 

In ref. |ji5|| , a particular ansatz B[D] was proposed to describe the first singularities 
close to the edge of the domain of convergence as a sum of "poles" with fractional powers 
(actually cuts). The free parameters were determined from the first orders of the perturba- 
tive expansion and the known properties of the Operator Product Expansion. This ansatz 
was used to compute D at arbitrary orders in perturbation theory, and then the integral 

using either FOPT or CIPT. In this case, when one increases the order, the evolution 
of the FOPT value of 5^ exhibits a plateau in agreement with the value obtained from 
the Borel sum D. 

This analysis is based on several assumptions. Because D and D share the same per- 
turbative expansion for real positive values of the coupling constant, the latter is expected 
to yield the " true" value of the Adler function for arbitrary (complex) values of the strong 
coupling constant in its convergence radius. One could then determine the "true" value of 
5° by integrating the Borel sum D(a(s)) over a circle in the complex s-plane, i.e. for com- 
plex values of a s . The ansatz for B[D] contains the three singularities that are relevant not 
only at high orders (where they are dominated by the first ultraviolet renormalons), but 
also at intermediate orders (where they are dominated by the first infrared renormalons), 
and even at low orders (since the first five orders of the perturbative expansion of D are 
used to determine the free parameters in B[D\). 
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In refs. [ 1Q| , 20, some aspects of this ansatz were discussed to map the Borel param- 
eter t into another parameter uu so that the cut plane along the real axis would be mapped 
into a disc of unit radius. It was argued that the series eq. fll.5| ) has optimal convergence 
properties once expressed in w, which would select this variable as the appropriate one 
to discuss renormalon models. Even though the main ingredients are the same (structure 
of the singularities, first orders of the perturbative series), the choice of w rather than 
t modifies the (non-singular) structure of the series eq. ( |1.5| ) and it was enough to alter 
significantly the outcome of the analysis. 

The present note aims at investigating other aspects of this ansatz. In Sec. 2, we recall 
the theoretical framework leading from the r decay width to the perturbative contribution 
5^°\ In Sec. 3, we present the ansatz used in ref. [15|, describing how the free parameters 



of the model are fixed. In Sec. 4, we discuss how the uncertainty on the Borel sum D can 
be underestimated when one use it to compute 5°, depending on the way the singularities 
related to infrared renormalons are avoided. In Sec. 5, we discuss in some detail the FOPT 
and CIPT prescriptions, and we identify two different ways of applying FOPT, leading to 
a further uncertainty for this prescription. In Sec. 6, we discuss how the free parameters of 
the model are fixed, based on a high-order expansion of the renormalon contributions which 
is affected by potentially significant corrections. We mimic the effect of such corrections by 
introducing a quadratic term in the series defining the Borel series, and study the impact 
on the FOPT/CIPT discussion. In Sec. 7, we study the notion of pole dominance, which 
was used in ref. fl5|] to favour their ansatz, and we see that several definitions could be 
imagined for such a dominance. In Sec. 8, we discuss the role of anomalous dimensions 
in the results of the discussion, in order to treat the first two infrared renormalons at the 
same level of detail. In Sec. 9, we show that the agreement of the FOPT or CIPT value 
with the Borel sum depends on the chosen weight, and is thus related to enhancement or 
cancellation of some parts of the integration contour. Finally, in Sec. 10, we summarise 
and conclude our study. 

2. Theoretical framework 

From the theoretical point of view, the r decay width R T can be described in terms of its 
contributions from non-strange vector ud, non-strange axial ud and strange us 

R T = R r y + R TA + Rr,s- (2.1) 
One can relate each of these to the corresponding spectral functions through 

-Ml / „ \ 2 



(2.2) 



R T = 12nS EW J T (l - (l + Im II« (s) + Im (s) 

and the spectral functions are related to the imaginary part of the correlators: 

= \v ud \^ v J\s) + n* J) (*)] + |K,| 2 [nri J )(.) + n^)(,)] (2.3) 
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combining two-point correlators of hadronic current with given angular momentum. They 
are defined from the correlator defined in QCD: 

KMP) = * J dxe^(0\TjX uD (x)jX uD (0)^\0) (2.4) 

where X = V or A, D = d or s, and the hadronic currents are J^ uD = D^/^u and 
jfiuD = D^itTi'oU. This correlator can be decomposed according to angular momentum: 

Since the correlator has only a singularity along the positive real axis, one can deform 
the integration contour and reexpress this quantity as 

K = 6„<w [(i + 2^) n« (s , + n<») (s) ] (2.6) 

= -mS EW i — (1 - ^) 3 kl + x)D^(M^x) + 4Z?(°)(M T 2 x)l (2.7) 
J\x\=\ x L J 

where we have rewritten the integral in terms of the Adler functions 

= -s^^(s)] D<°> = ^^[^(°)( S )] (2.8) 

If one performs the operator product expansion of D, one can write down an expression 
of R T of the form 

Rr,v/A = ^S EW \V ud \ 2 [l + «(°) 6 %a\ ( 2 -9) 



In the chiral limit considered to compute the perturbative contribution, = and one 
can consider only d( 1+0 ) which is independent of the renormalisation scale and has an 
expansion of the type: 



oo 



4 1+0) W 2 ) = ^E^Q (2-10) 

n=0 

with oq = a s (Q)/w. The first few orders are known: 

co,i = ci,i = i C2,i = 1.640 C3,i = 6.371 c 4 ,i = 49.076 (2.11) 

The value of cs ; i is still unknown, and a very frequent assumption consists in taking a 
geometric progression for these numbers, leading to C5 ; i ~ 283. 

The issue now consists in determining how to perform the integration over the contour 
for the perturbative contribution in the most accurate way. This issue is of importance 
here for two different reasons: 

• The main outcome of the analysis consists in the determination of a s (m 2 ), which 
is involved mainly through the perturbative contribution to the computation of the 
r-width 

• The Operator Product Expansion is expected to break down once one gets closer 
to the cut corresponding to the physical region, and thus its integration around a 
circular contour in the complex plane should be considered with a particular care 
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3. Model of renormalons 



In ref. [15], a renormalon model was presented to describe higher orders for the Adler 
function D eq. ( |2.10| ), and thus to compare different integration methods once higher order 
are taken into account. The ansatz consisted in one ultraviolet renormalon, corresponding 
to the sign-alternating divergence expected at higher orders in the perturbative series of 
D, and two infrared renormalons, mirroring the presence of condensates of dimension 4 
(gluon condensate) and 6 (higher order quark and gluon condensates) in the OPE of the 
Adler function taken in the chiral limit fl3|, [l6|, |l7|, [l8| . 

In the Borel plane, the ansatz for the Adler function has the following representation 

B[D](u) = B[DY v ](u) + B[Dl R ](u) + B[Di R }(u) + d PO + d PO u (3.1) 
where each renormalon "pole" has a cut singularity of the form 

d x 

B[D*}{u) = ( ^ )1+7 [1 + hip =F u) + b 2 (p T uf + . . .] (3.2) 

where the negative sign corresponds to an infrared (IR) renormalon, and the positive one 
to an ultraviolet (UV) renormalon. 

The corresponding model for the Adler function is given by 

poo 

D(a) = / dte-*/ a 5[D](t(u)) (3.3) 



with t = ttu//3q. The real part of this integral is expected to yield the "true" value of the 
perturbative series, whereas its imaginary part (divided by tt) should provide an estimate of 
the uncertainty, attached to the way one treats the singularities related to IR renormalons. 
The corresponding perturbative series can be worked out as outlined in sec. 5 of ref. [15] 

t?™ = pl+ 7rf +l0 f>+'+^(4)V (3 - 4) 

l p (6i + ci) / p\ 2 (b 2 + b 1 c 1 +c 2 ) fl_ 
P Q n + 7 \f3 J (n + 7)(n + 7-l) \n 3 

where ag = a s (Q)/-K. The values of the coefficients can be determined in the case of 
IR renormalons because their structure is connected with the contributions of the non- 
perturbative condensates occurring in the Operator Product Expansion of the correlator: 
p is related to the naive dimension of the condensate, 7 to its anomalous dimension, the 
b coefficients to the running of the strong coupling constant and the c coefficients to the 
perturbative series multiplying the condensate in the OPE. An analytic continuation can 
then be performed to assume a similar connection between UV renormalons and higher- 
dimension operators. 

The relative weight of the three renormalon contributions (indicated by dJ( v , , c^) 
was fixed using the 3rd, 4th and 5th orders in the expansion of the Adler function (03 1, C4 1, C5 1, 
the latter being set assuming a geometrical growth of the coefficients). All the coefficients 
c in eq. (|3.4j), related to the Wilson coefficients in the OPE, were set to zero apart from 
the coefficient c\ for the d 2 R pole. The first two terms (01,1,02,1) were not considered and 
are reproduced by adding an ad hoc term in the model (constant and linear terms in u). 
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4. Uncertainty on the Borel integral 



One can use eq. ( |3.3| ) in order to derive a "resummed value" of the perturbative expansion. 
The presence of IR poles on the positive real axis means that we have to give a prescription 
for the integral, depending on whether we take the contour integral above or below the real 
axis, which yields: 

b( s ) = b PV {s)±ib pole {s) (4.i) 

corresponding to the principal value and the pole contribution of the integral in eq. (|3.3|). 
When s is along the real axis, both Dpv and D po [ e are real, but this does not remain the 
case in the complex plane. 

We perform the contour integral in the complex plane eq. ( p.6| ) in order to compute 

5(0). 

<5 (0) = / ds K{s) D(s) = I ds K(s) [b PV (s)±ib pole {s)} (4.2) 

J\s\=M? J\s\=M% 

where Re Dpy (sq exp[i0]), Re D po i e (sQ exp[i(j)]) and Re K(sq exp[i^>]) are even functions 
of 4>, whereas the imaginary parts are odd functions of the same variable. Therefore the 
value of is obtained from 

= I ds K(s) Dp V (s) = I ds [Re K Re D PV — Im K Im D PV }(s) (4.3) 

J\s\=M? J\s\=M? 

The imaginary part of the integral evaluated with the principal value prescription (and 
divided by it) is sometimes taken as an estimate of the uncertainty on the value of the 



Borel integral [15]: 



Err 5^ = — (f ds K(s) D po i e ( s ) = <£ ds [Re K Re D po i e — Im K Im £> po i e ](s) 
K J\s\=M% J\b\=M* 

(4.4) 

Once several IR poles are included in the model, one must decide how to combine 
the contributions from D 1 ^, and in Err S^°K Ref. |ll| takes the sum of the two 

contributions with a relative sign maximising the error, i.e. 



Err = - 



7T 



/ dsK(s) [b I p 2(s)±bj^ e { 8 )] 
J\s\=M? 



(4.5) 



This amounts to assuming correlations for the variation of the uncertainty when one moves 
along the circle in the complex plane (the relative sign is assumed to be the same for any 
position along the circle). It seems more conservative to assume an absence of correlations, 
taking: 

Err 5(0) = - t ds K{s) [|^te(*)| + [D^MW (4-6) 
This prescription tends to increase the error bar on the Borel integral in a significant way. 
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5. Different methods of treating perturbation theory 



In refs. [15, [|, [7|, ^, 11] were discussed two main prescriptions to compute the integral 



eq. ( |1.3[) in terms of the perturbative expansion of the Adler function eq. ( 2.10| ), called 



Fixed-Order Perturbation Theory (FOPT) and Contour-Improved Perturbation Theory 
(CIPT). The methods were actually presented differently in refs.[||, 0, [8| and ref.[[l5|], and 
we will recall the salient elements of both presentations. 

The starting point of refs.Q, [8|, |j is the solution of the RGE which reads - as defined 



in [22, 



The full expressions for an arbitrary number of quark flavours (rif) are: (3q = \ (ll — f n /)> 

a 1 /ino 38„ \ a _ 1 / 2857 5033 ^ _i_ 325„2 N \ Qr „q 

Pi - 16 \ 102 ~ T n f) ' - 64 ^—2 IS -71 / + "ST™/ J ' and 



A = 266 



149753 _ / 1078361 6508 \ / 50065 6472 \ , 1093 , 

h 3564 C3 1 (3 )nt + 1 C3 H re? 

6 V 162 27 S J 1 \ 162 81 S V y 729 J 



where Ci={3,4,5} = {1-2020569, 7r 4 /90, 1.0369278} are the Riemann ("-functions. The stan- 
dard perturbative method to compute the contour integral consists then of expanding all 
the quantities up to a given power of a s (so). The starting point is the solution of the 
renormalisation group equation (RGE) for a s (s), which is expanded in a Taylor series of 
7] = ln(s/so) around the reference scale so J§ 

a s (s) = a s - f3 oV a 2 s + (-ftjj + /3 V ) a 3 s + ( - f3 2 rj + |/3 /3ir/ 2 - $jr/ 3 ^ a 4 

+ (-fin + -ftV + 3/3 /W - y /3 2 /3^ 3 + /3 V ) a* (5.2) 

+ ( - M + Ift/W + ^A)/V - yA^V ~ 6/3 2 /V + ^/3 3 /V " /3 V ) a° + O^ 6 ; oj) . 

Here the series has been reordered in powers of a s = a s (so) and we use the RGE j3- 
function ( |5.1|) . We insert eq.( |5.2|) up to a 5 s (since we know the (3 coefficients only up to 
that order) into the integral eq. ( |1.3|) , and we order the contributions according to their 
powers in a s , so that we obtain the familiar expression for fixed-order perturbation theory 
(FOPT) 0| 

00 

« (0) =E[««(0+5«(0]o?(^o), (5-3) 

n=l 

where the g n are functions of K m<n and /3 m < n -i, and of elementary integrals with loga- 
rithms of power m < n in the integrand. 

One can also perform the integral by performing a step-by-step integration along the 
circular contour, using the RGE solution eq. ( jsT^ ) to determine the value of the strong 
coupling constant at each point of integration. The initial point lies on the real axis (like 
for FOPT), and as one moves along the circular contour, the value of a s at a given point 
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is computed using RGE with the value obtained at the previous step. This second method 
was introduced in refs. (jj], |j) and is called Contour-Improved Perturbation Theory (or 
CIPT). Its advantages were discussed in ref. [11|. 

In ref. [JDJ, a different angle was chosen to present the methods. The starting point 
consisted not in inserting the solution of the RGE Q5.2|) into the perturbative expansion 



of the Adler function eq. (2.10), but rather in exploiting the fact that the Adler function 
is independent of the renormalisation scale. Indeed if we start from the perturbative 
expression for II: 



oo n+1 



n^s) = £ < £ c n, k f V = log ^ (5.4) 

^ — n i — n ' 



n=0 k=0 

with a M = a s (/i)/7r, we get the following expansion for the Adler function 



AT v nJr ^ 
n=0 fc=l 



In the language of ref. ]7], [| , we obtain for the function arising in eq. ( |5.3| 



ra+l 

fc=l 

But since the Adler function is independent of the renormalisation scale, as implemented 
in eq. ( |2.10| ), we can take the derivative of eq. (|5^ ) with respect to log// and exploit the 
RGE to reexpress the derivative of in terms of itself. We obtain an expansion in 
powers of and r] which is identical to zero. This yields expressions for c n ^>2 from c ny \. 
Since the values of /3 n >4 and c n ->§^ are unknown, the reconstruction of of c n ^>2 based on 
the RGE is only partial: the coefficients c n ^ with n > 6 and 1 < k < n — 4 cannot be 
computed fully because their equations involve some of these unknown coefficients. 

One can then perform the computation of the perturbative contribution to the tau 
width: 

« (,1 = EE fe »4f 1)3(1 +I)l0g *" 1 ( z ?)< (5J) 

n=lfc=i Lm J \ x \ = l X \ V / 

in two different ways, either by performing FOPT or CIPT. In FOPT, we set \i = M T , 
leading to 

oo n 

6 FOPT = E E kc n,k-h-l (5-8) 

n=l fe=l 

with 

1 f rlr 

Jk = ^-f — (l-x) 3 (l + x)log fc (-x) (5.9) 
In CIPT, we take fj 2 = -M^x to get 

4V = E c «.^n (5-10) 
n=l 
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with 



J* = J-I —(l -xf{l + x)a n {-M*x) (5.11) 



2m J\ x \ =1 x 



Since the two presentations are not obviously identical, it is interesting to determine 
if they are fully equivalent. In the case of CIPT, only the coefficients c n ^ coefficients from 
eq.( 2.10| ) are involved, so that the definitions are easily recognised as identical. In the case 



of FOPT, the situation is rather different, since we use the RGE at all orders to derive the 
coefficients c n ^>2 according to ref. fig] , whereas it is only exploited to up to a given order 
in the other references (a^ in ref. [Q]). The two results differ from each other through the 
higher-order terms that are kept or not. Indeed, the FOPT expression contains terms up 
to a certain order in a s {M T ) coming from: 

• the perturbative expansion of D, known up to an arbitrary order through the renor- 
malon model 

• the RGE of a s , used to reexpress the integral in terms of a series in a s (M T ), with 
coefficients known only up to /?3. 

The two previous definitions of FOPT can be rephrased in the following way: 



In ref. pq| , the RGE is used to expand a™(s) for n arbitrary, in powers of ^(sq) 
up to the required order in perturbation theory, setting all /3^>4 = 0. Then these 
expansions are put in eq. ( |2.10| ) to compute the FOPT value of 5(°> at a given order 
in perturbation theory. This method is denoted FOPT (B J) in the following. 



• In refs. Q, the expansion of a s , eq. Q5,2| ) up to O(af) is plugged into the series for 
5^°\ eq. (|T§|), the latter being cut at the required order in perturbation theory. We 
denote FOPT this way of dealing with higher orders. 

The difference between the two methods comes from a different use of the RGE for orders 
above the known terms (i.e., once /? n >4 is involved). For instance, for the n = 6 contribution 
to S^°\ we should use eq. ( |5.2[ ) up to a®, involving 04. This difference was already discussed 
in ref. pT| , where only low orders were known and included. In this reference, we discussed 
a method called FOPT + , which corresponds to the prescription denoted FOPT here, up to 
high-order coefficients in the perturbative expansion of D (modeled here, but set to zero in 
FOPT + ). A second method, called FOPT ++ in that paper, is related to the prescription 



denoted FOPT(BJ) here (in ref. |Tl|], we stopped the Taylor expansion at r\ 



n<5\ 



Following the methods presented in ref. [15|, we obtain the plot in Fig. |l] (in this paper, 
we will always use the illustrative value of the strong coupling constant a s (M T ) = 0.34). 
Let us emphasize that neither FOPT nor FOPT (B J) is "the true" FOPT as soon as we 
look for contributions from orders higher than n = 5, and they are both different from 
the complete value that would be derived by applying FOPT if the full RGE series were 
known: a piece is missing. 

One can see that the exploitation of the RGE at higher orders has a significant impact 
on 5(°\ The difference between the two methods remains within our (conservative) error 
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bars, but which exceeds the limited error bars chosen by ref. |15|| . In the case of CIPT, 
one has also to decide up to which order of perturbation theory one should write eq. (fT^), 
However, the integration of RGE is performed by small-step integration, and thus it is 
much less dependent on the exact cut on the power of r\ placed on eq. (|5.2|) 
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Figure 1: Value of 5^ as a function of the order and type of perturbation theory, corresponding to 
Fig. 7 from ref. fL5| . FOPT and FOPT(BJ) indicate two different ways of dealing with Fixed-Order 
Perturbation Theory as one goes to higher orders in perturbation theory. 



6. Extension of the renormalon model 

Based on the renormalon model presented in sec. [3|, ref. [15] fixed the renormalon residues 



di in eq. (3.1) by requiring the model to reproduce the first orders of the perturbative 
series 03^,04,1,05,1 (the latter being set to 05,1 = 283 following an ansatz concerning the 
geometrical growth of this coefficients). The argument stated in ref. |l5| is that only a 
constant term cLq° is actually needed to reproduce ci t i and 02,1 (even though a linear term 
df° is included, but turned out to be very small in the ansatz of ref. |T5[| ). 



With this assumption, ref. [15| shows that the balance between the two IR renormalons 



has a direct consequence on the discussion of CIPT/FOPT in their model: 

• A dominance of the d = 6 renormalon favours CIPT (in the sense that it yields a 
value close to the Borel sum of the series) 
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Pole 1 2 3 456789 10 11 12 



UV 36.1 11.7 3.7 1.6 0.8 0.5 0.3 0.2 0.1 0.1 0.1 0.1 
IR2 -377.4 46.2 7.5 2.5 1.2 0.6 0.4 0.2 0.2 0.1 0.1 0.1 
IR3 -106.3 -59.3 26.3 4.6 1.7 0.8 0.5 0.3 0.2 0.1 0.1 0.1 



Table 1: Relative contribution (in percent) Ep> / D™ x from 1/n 3, to a given order of perturbation 
theory for each pole of the renormalon model. The correction cannot be estimated for n = 0. 



• A dominance of the d = 4 renormalon favours FOPT (in the sense that it yields a 
value close to the Borel sum of the series) 

However, this procedure to fix the residues d, assumes that the perturbative expansion 
eq. Q3.4D is valid exactly when one neglects the contributions coming from the coefficient C2 
and from 1/n 3 remainders (since they are set to in ref. ) . Actually, these contributions 
can be quite significant, compared to other contributions that are explictly included in the 
renormalon model in ref [15]. Indeed, Tables [l] and |2| collect the relative contribution from 
1/n 3 and C2 terms: 



nd* 



p i+7r(i + 7 ) 

irdp 
p 1 +TT(l + 7) 



r(n + i + 7 ) 



r(n + 1 + 7) 



P 

±* 

P 



,n+l 



1 



n" 



(n + 7 )(n + 7 - 1) 



(6.1) 
(6.2) 



compared to the total contribution of a pole for a given order in perturbation theory as 
computed in ref IE]: 



D n p > X (a Q ) 



ird: 



x 



p 1 +rr(l+ 7 ) 



r(n + l + 7 ) ± 



A) 



,n+l 



(6.3) 



1 P (h+ci) 
f3 n + 7 



(62 + hci) 



(n + 7)(n -(- 7 - 1) 



In other words, the actual contribution to the perturbative expansion of the Adler series 
from a given pole at a given order would be 



,n+l[nn,X 

1 q y u P 



zE^ x + c 2 F^ x ] 



(6.4) 



z and C2 are unknown coefficients, in principle of order 1. One could neglect their presence 



at a given order a™ +1 of perturbation theory if Ep' /Dp' and Fp' /Dp' are small 
numbers (this is in particular assumed in ref [jn| for n = 2, 3,4). 

One notices that the perturbative coefficient for n < 3 is significantly affected by both 
kinds of contribution (E and F), and that the /-R3 pole is more affected by such corrections 
than /i?2- In addition, the value of C2 can affect significantly the situation for IR3 even for 



higher orders. Therefore, the model discussed in ref. [15] is likely to have its perturbative 
expansion affected by significant corrections at low orders - where the term "low orders" 
includes the 0(a 3 ) term (i.e., n = 2), which is used to determine the residues of the poles. 
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Pole 1 2 3 456789 10 11 12 



UV 36.4 28.1 6.7 2.8 1.5 0.9 0.6 0.4 0.3 0.3 0.2 0.2 0.1 
IR2 -29.7 -73.1 31.6 9.8 5.0 3.1 2.1 1.6 1.2 0.9 0.8 0.6 0.5 
IR3 -17.9 -23.7 -57.3 53.8 15.2 8.0 5.1 3.6 2.7 2.1 1.7 1.4 1.2 



Table 2: Relative contribution (in percent) F™ ,x / D™' from C2 to a given order of perturbation 
theory for each pole of the renormalon mode. 



Table 3: Coefficients of the ansatz for different values of d PO . 





d PO 


4° 


d% V 


4* 


4 R 


-9.116 


-3.834 


-1 


0.028 


7.82 


-176.8 


-4.167 


-1.913 


-0.5 


0.006 


5.49 


-95.16 


0.781 


0.008 





-0.0160 


3.16 


-13.53 


3.255 


0.968 


0.25 


-0.026 


1.99 


27.29 


5.729 


1.929 


0.5 


-0.037 


0.83 


68.10 


10.68 


3.850 


1 


-0.059 


-1.50 


149.7 



One can take into account the possibility of such large corrections at the lower orders 
by adding a quadratic term in u to the polynomial part, which will contribute to the 
perturbative expansion of the Adler function at O(a^) : 

B[D](u) = B[D^ v ]{u) + B[L>l R }(u) + B[b{ R ]{u) + d PO + df°u + d PO u 2 (6.5) 

The choice of d PO can be seen converted into a guess on the value of c^i, with the 
corresponding equivalence: 

d PO -1 -0.5 0.25 0.5 1 
~c 6 1291 2283 3275 3771 4267 5259 

The dependence of all results on d PO is linear, and one can perform exactly the same 
analysis as before, with the corresponding plots in Fig. As expected, O(l) values of d PO 
are enough to change the balance between the IR renormalon poles, as shown in Table [| 

One recovers the model in ref. |^] in the case where d PO = 0. The discussion of CIPT 
vs FOPT can be converted into a discussion on the value of d PO . The cases where d^ p 
vanishes (for d PO approximately between 0.5 and 1) correspond to cases where CIPT is 
preferred to FOPT if one wants an agreement with the value of 5^ from the Borel resum- 
mation. The cases where d\ R vanishes (approximately d PO between and 0.5) correspond 
to cases where FOPT is preferred to CIPT. 
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Figure 2: Value of 8^ as a function of the order and type of perturbation theory, corresponding to 
Fig. 7 from ref. FOPT and FOPTfBJ) indicate two different ways of dealing with Fixed-Order 
Perturbation Theory as one goes to higher orders in perturbation theory. The figures correspond 
to different values of the quadratic coefficient d^ = — 1, —0.5, 0, 0.25, 0.5, 1 (from left to right and 
top to bottom) 
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Table 4: Relative contribution (in %) to the coefficients of the perturbative expansion of the Adler 
function for different values of (if = —0.5,0,0.25,0.5 





Pole 


C4 1 






C7 1 


c 8 1 


Cg 1 


C10 1 


Cll 1 


C\2 1 


C13 1 




TJVi 


-3.8 


6.1 


-8.9 


9.2 


-15.9 


16.7 


-45.5 


33.7 


-576.4 


59.4 


-0.5 


IR9 


174.4 


236.1 


243.0 


164.8 


177.7 


112.6 


180.1 


77.1 


752.1 


43.7 




IR3 


-70.6 


-142.2 


-134.1 


-74.1 


-61.7 


-29.4 


-34.6 


-10.8 


-75.7 


-3.1 




UVi 


9.7 


-15.6 


15.8 


-38.6 


30.3 


-236.7 


54.0 


200.0 


77.6 


119.5 





IR 2 


100.4 


135.9 


97.5 


155.9 


76.3 


360 


48.3 


-103.6 


22.9 


-19.9 




IR 3 


-10.0 


-20.2 


-13.3 


-17.3 


-6.5 


-23.2 


-2.3 


3.6 


-0.6 


0.3 




UVi 


16.4 


-26.5 


23.3 


-96.5 


45.6 


314.6 


72.3 


127.2 


89.4 


107.0 


0.25 


IR 2 


63.3 


85.7 


53.4 


145.0 


42.7 


-178 


24.1 


-24.5 


9.8 


-6.6 




IRs 


20.2 


40.8 


23.3 


51.4 


11.7 


-36.7 


3.6 


-2.7 


0.8 


-0.4 




UVi 


3.1 


-37.4 


29.0 


-258.7 


57.8 


159.4 


84.2 


110.4 


95.4 


102.5 


0.5 


IR 2 


26.3 


35.6 


19.6 


114.6 


15.9 


-26.6 


8.3 


-6.3 


3.1 


-1.9 




IR 3 


50.5 


101.8 


51.3 


244.1 


26.3 


-32.9 


7.5 


-4.1 


1.5 


-0.6 



Larger (positive or negative) values of d^° correspond to values where neither FOPT 
nor CIPT yield a plateau with a value in agreement with the Borel resummation. CIPT 
yields a result that is stable, but in disagreement with the value obtained by Borel resum- 
mation, whereas FOPT yields an unstable result which sometimes crosses the Borel result. 
We see that the value of the Borel resummation depends significantly on the value of cZf" > 
as well as the uncertainty attached to it. 



7. Definition of pole dominance 

From the previous section, one can see that a seemingly small change in df° has an 
important impact on higher orders and on the value obtained from the Borel resummation. 
It seems difficult to determine a priori which model is more relevant, and different criteria 
can be imagined. A rather usual approach consists in assuming that the first IR pole should 
"dominate" over the following ones for the model to be reasonable. Let us remark first that 
such a requirement is by no means mandatory in the framework of the renormalon approach, 
in which the relative contributions form different poles is free, at least in principle. For the 
time being, let us assume that such a dominance is indeed the sign of a good model. The 
next question is: what is the definition of this dominance in practice ? 

A first definition, chosen in ref. [{15}j, consists in considering the contribution of the 
different poles to the expansion of the Adler function at intermediate orders (for n be- 
tween 4 and 8). A dominant pole gives the largest contribution to each coefficient of the 



perturbative series for intermediate orders. In ref. [15], the case d^° = is considered 
as interesting because the relative contribution to a given c n , 1 from IR2 is larger than 
IR3. This is recalled in Table |J Let us remark that already for n = 7, large cancellations 
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Table 5: Order-by-order results for 8^ FOPT(BJ) (top) and CIPT (bottom), for d%° = -0.5 



Order 


1 


2 


3 


4 


5 


6 


UVi 


-0.00682 


-0.01114 


-0.01277 


-0.01281 


-0.01219 


-0.01161 


IR 2 


-0.90568 


-1.28083 


-1.37745 


-1.33435 


-1.25036 


-1.17583 


ID 

IK3 


2.467 iu 


3.5828b 


no A TO 

3.99478 


3.99000 


3.82183 


3.b4055 


Pol 


-1.41698 


-2.12173 


-2.40198 


-2.42283 


-2.33050 


-2.22195 


Sum FOPTfRTl 


10822 


16916 


20258 

I J . Zj VJ Zj O O 


22001 

\J . Z. Zj 1 


22878 


231 1 6 

U. ZO 1 1U 


UVi 


-0.00932 


-0.01194 


-0.01141 


-0.01174 


-0.01150 


-0.01165 


IR 2 


-1.23768 


-1.27781 


-1.26035 


-1.24537 


-1.23643 


-1.23238 


IR 3 


3.33131 


3.64857 


3.67319 


3.66712 


3.66174 


3.65950 


Pol 


-1.93642 


-2.18124 


-2.21166 


-2.21166 


-2.21166 


-2.21166 


Sum (CIPT) 


0.14790 


0.17758 


0.18977 


0.19836 


0.20214 


0.20381 


Borel sum 


0.26036 













Table 6: Order-by-ordcr results for 5^ FOPT(BJ) (top) and CIPT (bottom), for d%° = 0. 



Order 


1 


2 


3 


4 


5 


6 


UVi 




0.01734 


0.02836 


0.03248 


0.032583 


0.031007 


0.029536 


IR 2 




-0.52120 


-0.73709 


-0.79269 


-0.767886 


-0.719556 


-0.676666 


IRs 




0.34660 


0.50942 


0.56798 


0.567304 


0.543394 


0.517619 


Pol 




0.26548 


0.36847 


0.39480 


0.388012 


0.373940 


0.362262 


Sum 


FOPT(BJ) 


0.10822 


0.16916 


0.20258 


0.22001 


0.22878 


0.23275 


UVi 




0.02370 


0.03039 


0.02903 


0.02986 


0.02927 


0.02965 


IR 2 




-0.71226 


-0.73534 


-0.72530 


-0.71668 


-0.71154 


-0.70920 


IRs 




0.47365 


0.51876 


0.52226 


0.52140 


0.52063 


0.52031 


Pol 




0.36280 


0.36378 


0.36378 


0.36378 


0.36378 


0.36378 


Sum 


(CIPT) 


0.14790 


0.17758 


0.18977 


0.19836 


0.20214 


0.20453 



Borel sum 0.23709 



occur between the UV pole and the first IR pole. If one now considers different values of 
d,2°, one notices that d^ = 0.25 is also in good agreement with this criterion. In this 
case, there is no plateau for FOPT at intermediate orders in perturbation theory, and the 
outcome of the Borel resummation lies between FOPT and CIPT. Larger positive values 
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Table 7: Order-by-order results for 5^ FOPT(BJ) (top) and CIPT (bottom), for d%° = 0.25 



Order 


1 


2 


3 


4 


5 


6 


UVi 


0.02943 


0.04811 


0.05511 


0.05528 


0.05260 


0.05011 


IR 2 


-0.32896 


-0.46522 


-0.50031 


-0.48466 


-0.45415 


-0.42708 


J±t3 


-u.oyoyo 


-l.UZ 1 oU 


-1.14041 


i i/i/iri/i 

-1.144U4 


-i.uyoo/ 


-1.U4O50 


Pol 


1.10671 


1.61357 


1.78311 


1.77170 


1.70077 


1.63356 


Sum 


0.10822 


0.16916 


0.20258 


0.22001 


0.22878 


0.23355 


UVi 


0.04021 


0.05156 


0.04925 


0.05066 


0.04966 


0.05030 


IR 2 


-0.44954 


-0.46412 


-0.45778 


-0.45234 


-0.44909 


-0.44762 


IRs 


-0.95518 


-1.04615 


-1.05320 


-1.05146 


-1.04992 


-1.04928 


Pol 


1.51241 


1.63629 


1.65150 


1.65150 


1.65150 


1.65150 


Sum 


0.147898 


0.17758 


0.189766 


0.198358 


0.202144 


0.204897 


Borel sum 


0.22546 













Table 8: Order-by-order results for FOPT(BJ) (top) and CIPT (bottom), for df° = 0.5 



Order 


1 


2 


3 


4 


5 


6 


UVi 


0.04151 


0.06787 


0.07774 


0.07797 


0.07420 


0.07068 


IR 2 


-0.13672 


-0.19335 


-0.20793 


-0.20143 


-0.18875 


-0.17750 


IRs 


-1.74451 


-2.56403 


-2.85881 


-2.85539 


-2.73504 


-2.60531 


Pol 


1.94795 


2.85867 


3.17143 


3.15538 


3.02761 


2.90485 


Sum FOPT(BJ) 


0.10822 


0.16916 


0.20258 


0.22001 


0.22878 


0.23435 


UVi 


0.05672 


0.07273 


0.06947 


0.07146 


0.07004 


0.07095 


IR 2 


-0.18683 


-0.19289 


-0.19026 


-0.18799 


-0.18664 


-0.18603 


IRs 


-2.38401 


-2.61105 


-2.62867 


-2.62433 


-2.62047 


-2.61887 


Pol 


2.66202 


2.90880 


2.93922 


2.93922 


2.93922 


2.93922 


Sum CIPT 


0.14790 


0.17758 


0.18977 


0.19836 


0.20214 


0.20526 


Borcl sum 


0.21382 













of d,2° would be disfavoured according to this definition of dominance. On the other hand, 
negative values of in particular large ones, fulfill this definition, as can be seen for 
instance from Table. |j. 

But one may wonder whether the relative contribution to intermediate orders of D 
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the best way of determining whether one element of the model or another yields the most 
significant contribution to S^ '. Indeed, the intermediate orders contribute only mildly to 
the actual value of 5^°', since they are multiplied by higher and higher powers of a s . It 
seems rather natural to break down the contribution to 5^ into the contributions from 
the UV pole, the two IR poles and the polynomial term. 

An alternative definition of the dominance from one pole would correspond to providing 
most of the contribution to 5^ . Looking at Table 0, we see that in the case 
0, and contrary to our intuition, the so-called dominant IR2 pole yields a contribution that 
is almost canceled by the so-called subdominant IR3 pole - the contributions from the 
two IR poles and the polynomial part being of the same order of magnitude. When one 
increases (see Tables |7| and ^) , one can notice that the contribution from the IR3 pole 
grows and is canceled by the polynomial part. The value d^° — 0.08 corresponds to the 
extreme situation where the residue of the IR3 pole vanishes (and thus this pole does not 
contribute) and the IR2 pole saturates the contribution from IR renormalons. It seems fair 
to require all the various contributions (individual "poles" and polynomial term) to yield 
contributions of the same size, which is the case for -0.5 < dl° < 0.25 (see Tables 
corresponding to a rather wide range of behaviours of FOPT/CIPT/Borel sum. 

8. Anomalous dimension for the operator of dimension 6 



As can be seen from eq. (|3.2j), the term of renormalon "pole" is slightly misleading, since 
one expects radiative corrections to turn these poles into cuts in the Borel plane. One 
can relate both types of renormalons to QCD operators. In particular, the presence of 
IR renormalons mirrors the existence of condensates in the OPE of the correlator under 



scrutiny [17]. In particular, the anomalous dimensions of the latter are used to fix some 



unknown parameters of the model. Ref. |15|| considers: 

• for UV\, a vanishing anomalous dimension for a single "effective" condensate 

• for IR2, the anomalous dimension corresponding to the gluon condensate [eq. (5.13) 
in this reference] 

• for IR3, a vanishing anomalous dimension for a single "effective" condensate 



As explained in ref. [17|, the structure of the cuts becomes rather involved once the full 
set of operators are considered. IR2 is rather simple since only one operator is involved 
by dimensional arguments: this renormalon is linked to d = 4 operators, and only the 
gluon condensate is involved (neither the identity operator nor the quark condensate since 
we work in the massless limit). On the other hand, both UV\ and IR3 are related to 
dimension-6 operators, namely: 

(V^V0(V>7 M V0 (V>7^75 VO (VV75VO (8-1) 
OH^)(^r A vo (^75^X^75^) IabcG^g; b g^ c (8.2) 



In ref. [17], the anomalous dimensions of d = 6 operators were reconsidered. After 



diagonalising the RG mixing matrix, the diagonal operators were shown to have anomalous 
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Figure 3: Figures for different values of the anomalous dimension of the effective d = 6 condensate: 
7 W = Pi\ with A = 0.379 (left) and -0.753 (right). 



dimensions at one loop of the form 

7 W = Ai = {0.379, 0.126, -0.332, -0.753, 0} (8.3) 

These five contributions should a priori be included individually in the renormalon model, 
rather than through a single "effective" condensate of unclear anomalous dimension, and 
actually set to zero. Since the d = 4 operator is described with its correct anomalous 
dimension and since we compare d = 4 and d = 6 renormalons, it seems fair to treat both 
renormalons on the same footing. 

As can be seen from the perturbative expansion eq. ( |3.4[) , the large-order behaviour of 
a "pole" is: 

4 1+0) (Q 2 ) = E r " a Q r n (^) n n\ (8.4) 

n=0 \ P / 

We see that the larger the anomalous dimension 7, the more dominant the pole at large 
n. But conversely, at smaller n, the operators with smaller anomalous dimensions compete 
(and can even be more significant) than the pole with the largest n. Therefore, once the 
proper cut structure of the second IR pole is taken into account, it becomes necessary but 
difficult to fix the relative "strengths" of the cuts from the different operators from the 
lowest orders of perturbation theory. 

As an illustration of the role of anomalous dimensions in the discussion, we set = 
PlA with A = —0.753 and 0.379 (rather than in the original model), and follow the same 



procedure as in ref. [15| to obtain the two plots in fig. y, indicating a rather wide range 
of behaviours for FOPT, depending on the choice of anomalous dimensions (CIPT on the 
other hand remains very stable). In addition, the resummed values for 6^ can vary from 
0.23455 to 0.25553 if one changes this single parameter. 
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Table 9: Relative contribution (in %) to the coefficients of the perturbative expansion of the Adler 
function for different values of ce,i, for a model with two IR poles for u = 3. 



Ca 1 


Pole 


C/i 1 






C7 1 

w ,1 


Cq 1 


Cq 1 






fin 1 




3275 


UVi 


9.7 


-15.7 


15.9 


-38.9 


30.4 


-241.3 


54.1 


199.1 


77.7 


119.4 






100. 


136.6 


98.0 


157.3 


76.7 


367.5 


48.5 


-103.3 


23.0 


-19.8 




IR3a 


-9.4 


-12.8 


-8.2 


-11.1 


-4.4 


-16.7 


-1.7 


2.8 


-0.5 


0.3 




1^36 


-1.2 


-8.0 


-5.6 


-7.3 


-2.7 


-9.4 


-0.9 


1.3 


-0.2 


0.1 


2283 


UVi 


5.8 


-9.4 


13.6 


-59.6 


56.9 


95.6 


139.3 


77.8 


128.0 


87.1 




IR 2 


-81.4 


-110.3 


-113.5 


-324.6 


-193.3 


-196.0 


-168.2 


54.3 


-50.9 


19.5 




IR 3 a 


179.4 


244.9 


225.2 


542.9 


262.9 


-211.3 


140.9 


-34.9 


24.7 


-7.1 




R^3b 


-3.8 


-25.3 


-25.4 


-58.7 


-26.4 


19.6 


-12.1 


2.8 


-1.8 


0.5 


4267 


UVi 


13.6 


-22.0 


17.0 


-33.8 


25.3 


-96.2 


42.9 


599.3 


66.5 


142.0 




IR 2 


283.1 


383.2 


211.0 


274.4 


128.2 


293.7 


77.1 


-623.0 


39.4 


-47.4 




IR 3 a 


-198.0 


-270.4 


-133.0 


-145.8 


-55.4 


-100.5 


-20.5 


127.0 


-6.1 


5.5 




R^3b 


1.4 


9.1 


4.9 


5.2 


1.8 


3.1 


0.6 


-3.3 


0.1 


-0.1 



We can extend the model of ref. [15] by assuming a value for c^i and splitting the 



second IR pole into two different poles, with different anomalous dimensions: 

B[b)(u) = B[DY V ](n) + B[b{ R ]{u) + B[D™](u) + B[d£](u) + d™ + d?°u (8.5) 

We denote IR^a for the cut with A a = 0.379 and IR^b with = —0.753. If we vary C6,i 
between 2283 and 4267 (values corresponding d% = ±0.5 in the previous model) and look 
at the relative contribution of each pole to a given order of D (Table ||), we see that for 
values above 3000, one has a dominance of IR2, with a significant cancellation of IR^ a an d 
a growing contribution from IR^b- On the other hand, below 3000, IR^ a takes over IR2, 
whereas IR3b remains small. If we now consider the contributions to 5^ (Tables 10|- 1^) 



and if we require them to be of similar order, values of C6,i above 3000 are acceptable. In 
this case, it is quite interesting to notice that the breakdown in terms of the contributions 
from different poles is very different, as well as the relative contributions to given orders of 
perturbation theory. Even though the value of 5^ is quite similar in fine, the agreement 
with FOPT and CIPT depends quite strongly on the value chosen for cgj. 

9. Other moments 

One can use the same machinery to analyse further moments, for instance some higher 
moments used to determine the condensates from a fit to r spectral data (fig. ||) J?], H (|, 11 
and moments proposed in ref. |2|| (fig.||): 

R t,v+a= 2 / ds 1 2 ( — J [ 1 + —)(M s ) + ai(s)) 
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Table 10: Results for the FOPT(BJ) integral, for c 6 ,i = 3275 in a model with two IR poles for 
u = 3. 



Order 


1 


2 


3 


4 


5 


6 


UVi 


0.01740 


0.02846 


0.03259 


0.03269 


0.03111 


0.02964 


IR 2 


-0.52346 


-0.74029 


-0.79613 


-0.77122 


-0.72268 


-0.67960 




0.03761 


0.05553 


0.06153 


0.06031 


0.05618 


0.05208 




0.41127 


0.59763 


0.66170 


0.65995 


0.63438 


0.60771 


Pol 


0.16540 


0.22783 


0.24289 


0.23829 


0.22979 


0.22293 


Sum FOPT(BJ) 


0.10822 


0.16916 


0.20258 


0.22001 


0.22878 


0.23275 



Borel sum 0.23655 



Table 11: Results for the FOPT(BJ) integral, for c 6 ,i = 2283 in a model with two IR poles for 
u = 3. 



Order 


1 


2 


3 


4 


5 


6 


UVi 


0.01044 


0.01707 


0.01955 


0.01961 


0.01866 


0.01778 


IR 2 


0.42298 


0.59818 


0.64330 


0.62317 


0.58395 


0.54915 




-0.72479 


-1.07018 


-1.18572 


-1.16230 


-1.08269 


-1.00363 




1.29039 


1.87510 


2.07611 


2.07060 


1.99040 


1.90671 


Pol 


-0.89078 


-1.25101 


-1.35066 


-1.33107 


-1.28154 


-1.23885 


Sum FOPT(BJ) 


0.10822 


0.16916 


0.20258 


0.220014 


0.22878 


0.23116 


Borel sum 


0.19598 













tin 6|V^| 2 Sew 
u t,v+a — ■ 



mi 



ds 



N s 1 
+ 



N-lml N-l\m 



N 



As in the case of the r width, we can use OPE to expand these moments. If we focus on 
the perturbative contribution (D = 0), we can reexpress it as an integral over the circle 
using integration by part: 



1 + 5 i0) = -2m <f —w(s)[D(s)] D=0 

J\s\=s s 



(9.1) 



2m I dsu(s)[U^} D=0 {s) = 4tt 2 f° u(s)[lm n( 1+0 )(s)] D=0 

^|s|=s Jo 



\s\=so 

where w(s) = ds' u(s'), so that we have 



u i(s) 



i so 

2 

so 



so 



1 + 2- 



soj \so 



(9.2) 
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Table 12: Results for the FOPT(BJ) integral, for c§ : \ = 4267 in a model with two IR poles for 
u = 3. 



Order 


1 


2 


3 


4 


5 


6 


UVi 


0.02437 


0.03984 


0.04563 


0.04577 


0.04356 


0.04149 


IR 2 


-1.46990 


-2.07876 


-2.23557 


-2.16562 


-2.02932 


-1.90836 




0.80001 


1.18125 


1.30878 


1.28293 


1.19506 


1.10779 




-0.46784 


-0.67983 


-0.75271 


-0.75071 


-0.72163 


-0.69129 


Pol 


1.22158 


1.70666 


1.83645 


1.80764 


1.74112 


1.68471 


Sum FOPT(BJ) 


0.10822 


0.16916 


0.20258 


0.22001 


0.22878 


0.23434 


Borel sum 


0.27712 













wqi{s) 



12 



(I + 1)(Z + 3)(Z + 4) I + 1 \s J I + 3 \s 



1+1 6 fs\ l+3 4 / S ^+ 4 



+ 



1 + 4 \s 



- (9-3) 



uu(s) 



wu(s) 



-1 i-JL 

SO V s 



1 + 2 " " 

so J \so 



12(3/ + 7) 



(Z + l)(Z + 2)(/ + 3)(/ + 4) l + l\s 

1+3 nri / _ \ i+4 



+ 



6 



1 + 3 



10 

Z + 4 



1+1 



+ 



+ 



1 + 5 



1 + 2 \s 
so 



1+2 



(9.4) 
(9.5) 



u N {s) 
w N {s) 



1 

so 



N 



N-ls 



+ 



1 



N s 
'2{N + l) + V 



N - 1 

N 

2(N - 1) 



AT' 



s 



1 



JV+l 



(9.6) 
(9.7) 



The moments utv were introduced in ref. |2J] to suppress the higher dimensional conden- 
sates that were noted to affect the analysis of the pinched weight moments in refs. 0, [8]. 

[TT[1 . These moments uat were used to extract the strong coupling constant by fitting 
the tau data, once the quark and gluon condensate were set to fixed values. The authors 
extracted the information on two different quantities (the strong coupling constant and a 
high-dimension condensate) by fitting the integrals I n (sq) obtained with the same weight 
un but different radii for the contour of integration so (between 2.3 GeV 2 and m 2 ). The au- 
thors claimed an impressive agreement between the values of the strong coupling constant 
obtained for different N. 

Such an agreement is not particularly surprising. Let us first of all notice that the 
points from r data between s* = 2.3 GeV 2 and m 2 . are correlated and have significant 
uncertainties, meaning that the input for the fit is essentially one integral, say I N (s*), 
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the integrals for other values of so carrying very little additional information. The fit is 
therefore perfect, with one input and two parameters (the strong coupling constant and a 
high-dimensional condensate). Moreover, the output of the fit is indeed very stable as far 
as the strong coupling constant is concernend, since these weights can be rewritten as: 



un(s) = — 
so 




(9.8) 



Once inserted in the integral used to computed 5^°\ and taking only power corrections 
(without logarithms) for the the OPE of II, one can see that the first bracket provides es- 
sentially a correlation between a s and the dimension-four condensates, and this correlation 
is identical for all the values of N. The second one fixes the value of the condensate of 
dimension 2N + 2 in terms of the dimension-4 condensate. Two sum rules for two different 
values of N provide therefore the same correlation between a s and the gluon condensate, 



fixed in the analysis of ref [24|. 

When we compare the figures [| and [| for the different moments, it is not clear whether 
CIPT or FOPT should be preferred in such a context. The moments tend to put a different 
emphasis between the contribution from d = 4 and d = 6 renormalons, which alter the 
discussion followed previously for 6^°'. It is quite interesting to notice that the agreement 
between FOPT and the Borel sum is not automatic, and depends on the structure of the 
kernel considered. In order to quantify this, one can consider the difference between the 
contribution of the IR poles in FOPT/CIPT and their Borel resummed values, as shown 



in Table 13 for Moo; ^io> an d u%. The agreement with the Borel resummed version for 
n around 7 is better for CIPT in the case of u%, equally bad for FOPT and CIPT in the 
case of Uki for k = 1, 1 = 0, and better for FOPT in the case of k = 1, 1 = 1. 
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Figure 4: Figures for Uki moments with fc = 1, 1 — (left) and k — 1, 1 — 1 (right). 
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Table 13: Results for the discrepancy between contributions for IR poles in FOPT/CIPT com- 
putations and the Borel sum for various weights 
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Figure 5: Figure for un moment with N = 2. 



10. Conclusion 

In this paper, we have investigated several aspects of renormalon models for the Adler 



function, used recently to compare fixed-order and contour-improved perturbation theories 
(FOPT and CIPT) to treat the contour integral for the theoretical estimate of the r 
width ^(J . Indeed the difference between the two treatments induce a significant 
systematics on the extraction of the strong coupling constant at the r mass. The particular 
renormalon ansatz for the Adler function in ref. [|l5|] suggested that FOPT was to be 
preferred with respect to CIPT, since it converges to the value of the Borel sum (taken as 
the true value of the integral). 

During our study, we have noticed the following points: 

• Once several infrared poles are included in the model, one needs to define how one 
combines the uncertainties estimated from the Borel integral. Depending on the 
treatment of the RGE for a s at higher orders (for which the f3 function is not known) , 
the FOPT prescription can yield noticeable differences at intermediate orders, in 
agreement with more conservative estimates of the uncertainty. 



• Ref. [15 1 sets to zero the contributions from unknown terms in the Wilson coefficients 
involved in the OPE of the Adler function . Moreover, the perturbative expansion 
of the ansatz is obtained through an expansion in powers of 1/n (where n is the 
order of perturbation theory) where only the first order are kept. This (truncated) 
perturbative expression is used to determine the residues of the poles from the first 
order of the perturbative expansion of D. We noticed that the neglected contributions 
are not particularly small at the orders of perturbation theory used to determine these 
residues, which can affected by significant uncertainties. 

• We extended the renormalon models by taking into account such potentially large 
contributions in particular for n = 3, and we investigated some cases where FOPT or 
CIPT are in better agreement with the values obtained from the Borel integral. We 
discussed two different definitions of the dominance of a pole, in order to determine 
which cases of these extended models could be considered as acceptable because of 
the dominance of the first infrared pole. 

• We examined the issue of the anomalous dimension of the second infrared pole, which 
actually corresponds to five operators of different dimensions, and observed rather 
different behaviours of the perturbative series. 

• We discussed other weights, noticing that the better agreement of FOPT with the 
value of the Borel integral is not a universal feature, and depends on which part of 
the contour integral is suppressed or enhanced by the weight. 

Renormalon models provide very attractive features to discuss qualitative aspects of 
higher-order perturbation theory. However, in the present discussion, we want to compare 
small differences between two treatments of perturbation theory, requiring a quantitative 
model of the higher orders of the Adler function. Any given ansatz based on renormalon 
calculus involve a large number of unknown coefficients both for the singular terms (residues 
of the renormalon "poles") and the non-singular terms (polynomial contribution). Only a 
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limited number of these coefficients can be fixed through the first few known terms of the 
perturbative expansion of the Adler function - the other ones being generally set to zero. 
It is not clear that the simplified description of the renormalon singularities by an ansatz, 
assumed to be valid at high orders, is sufficient at such low orders. 

The particular ansatz chosen in ref. [15] does not exhaust the potentialities of model 



building provided by renormalon calculus, and we have described a few extensions leading 
to rather varied conclusions concerning the comparison of FOPT versus CIPT. Our study 
shows that this particular ansatz cannot be taken as part of a reference test to determine 
whether FOPT, CIPT or yet another method should be adopted to extract a s from hadronic 
r decays. Moreover, significant systematics (as large as the difference between the standard 
CIPT and FOPT results) ought to be added to the results based on such an ansatz, since 
it is only one among many different renormalon models for perturbative expansions at high 
orders. 

The previous discussion is essentially based on the fact that we assume the Borel sum 
eq. ( |Q| ) to provide the "true" value of the asymptotic perturbative series for the Adler 
function. In particular, its value is used to determine whether FOPT or CIPT should be 
preferred. Let us mention that the theoretical estimation of the r decay is rather particular 
in this respect, since the low value of the r mass compared to hadronic scales requires one 
to compute an integral over a contour in the complex energy plane. The assumption that 
the Borel sum yields the true value of the Adler function should hold not only for real 
positive values of the coupling constant, but also for values of a s (s) where s is complex. 
The theory of asymptotic expansions |^q] indicates that there are functions for which this 
continuation of asymptotic series is not simple and one might encounter discontinuities 
when one crosses frontiers in the complex plane (Stokes lines). It would be interesting to 
determine whether such a situation could occur in renormalons models, and what their 
impact could be in the issues discussed here. 
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